Creating quantum correlations through local non-unitary memoryless channels 
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We show that two qubits, initially in a fully classical state, can develop significant quantum correlations as 
measured by the quantum discord (QD) under the action of a local memoryless noise (specifically we consider 
the case of a Markovian amplitude-damping channel). This is analytically proven after deriving in a compact 
form the QD for the class of separable states involved in such a process. We provide a picture in the Bloch sphere 
that unambiguously highlights the physical mechanism behind the effect regardless of the specific measure of 
QCs adopted. 
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The existence of states where two or more systems are cor- 
related in a way unattainable in classical physics ranks among 
the most puzzling and yet distinctive features of quantum me- 
chanics. Such possibility is commonly pictured in terms of 
an extra amount correlations, usually referred to as quantum 
correlations (QCs), which a multipartite system can possess in 
addition to those of a merely classical nature. Until recently, 
the scrutiny of QCs has been almost ubiquitously intertwined 
with investigations on entanglement and the pivotal role 
that it plays in the area of quantum information processing 
A breakthrough yet occurred as soon as it was realized [jsj 
that while classicality always entails separability the reverse is 
in general untrue (a state is entangled iff it is non-separable). 
Entanglement thereby is not the only form in which QCs can 
occur. Such finding brought about a widened perspective, 
which is currently prompting a growing number of researchers 
to advance the field along various lines. A prominent one is 
the quest for faithful easy-to-handle indicators of QCs flS]. 
Among those proposed so far, quantum discord (QD) yj is 
having a considerable impact despite its explicit calculation 
is usually demanding even for two qubits (i.e. a pair of two- 
dimensional systems). Yet, evidence of its ability to cap- 
ture QCs not detected by entanglement has been supplied in 
various frameworks such as one-qubit quantum computation 
(even experimentally) |01 and quantum phase transitions ( H ] . 

Another major concern that soon arose is to assess how 
QCs, according to such a novel paradigm, are affected by 
non-unitary dynamics. These typically stem from the inter- 
action with an environment, a process where entanglement 
is extremely fragile in most cases ill 01 ■ In contrast, QD 
was proven to be in general quite resilient to such dynamics 
and, strikingly, in some cases even fully insensitive over long 
stages ITj]. From a reverse perspective, it was shown that the 
quantum noise arising from a common bath can create QCs 
initially fully absent 18|], a phenomenon well-known to occur 
for entanglement either Also, preexisting QCs can exhibit 
an increasing behavior at some stages of their time evolution 
in the presence of non-Markovian local channels as recently 
demonstrated for both entanglement and QD Holl . In such in- 
stances, the increase of QCs stems either from the ability of a 
common memoryless reservoir to mediate an effective inter- 
action or as a memory effect of local environments. The ques- 
tion is now raised: Can the interaction with a bath which is 
both local and memoryless enhance QCs? The answer is well- 



known to be negative for entanglement, which cannot grow 
under any local quantum maps yj] . As for the full amount of 
QCs, however, the issue is not as much trivial since as stressed 
above even separable states may feature some quantumness 
As far as local unitaries are concerned, though, QCs' 
measures, such as QD, cannot increase IShS]. In some re- 
spects, this conclusion might be expected to be strengthened 
with noisy local operations. However, QCs without entan- 
glement can arise under a local non-unitary operation owing 
to its ability to map orthogonal into non-orthogonal states Isft. 
In this work, we present a simple paradigmatic process, where 
entanglement is absent throughout, which clearly testifies that 
QCs can even be entirely created (or increased) solely via the 
interaction with a bath which is both local and memoryless. 
We provide rigorous and comprehensive insight into such ef- 
fect in a way that makes transparent the underlying physical 
mechanism. Although we focus on QD, it will become clear 
that the essential physical effect takes place regardless of the 
specific measure of QCs used. 

To begin with, we briefly recall the definition of QD. Given 
two systems A and B in a state p, this measures the discrep- 
ancy between the mutual information I and the classical cor- 
relations C associated with p fsl. A local measurement on B in 
a given orthonormal basis can be specified by a complete set 
of projectors [Bk], where k indexes a possible outcome. If k 
is recorded with probability /?a =Tr[BApBi.] the overall system 
collapses onto the (normalized) state pt - {B^pBk)! Pk- Then 
the QD can be expressed as fsj] 



D^(p) ^S(pb)-S{p)+ min V pkS(pk) 



(1) 



In Eq. pb - TtaP is the reduced density operator that de- 
scribes the state of B, S (cr) - -Tr (p log20-) is the Von Neu- 
mann entropy of an arbitrary state cr and the infimum in the 
last term is evaluated over all the possible sets {Bk}- As for the 
QD involving measurements on A with associated projec- 
tors {Ak}, this is obtained from ^ through replacement of B 
with A. In general, D^^D^ 

Specifically, in the process under study A and B are two 
qubits and the only involved states have the separable form 



i(|0)A<0|®T0B- 



|1)a<1|®tib) 



(2) 



where to(i) are generic density matrices while the orthonor- 
mal set of vectors {|0), 1 1)) is the usual local computational ba- 




FIG. 1. (Color online) Behavior under a local amplitude-damping 
channel of discord D (a) and classical correlations C (b) against p. 
In each panel, the inset shows the corresponding behavior expressed 
as a function of the rescaled time yt, where p= 1-e"'" . 



sis. The states (|2]) are commonly dubbed as classical-quantum 
states owing to = as is immediate to check 1121], while 
in general i^Q. We will thereby set DsD^ henceforth. 
Consider now the initial state 



po = ^ (|0)a<0| ® \+}b{+\ + |1>A<1| ® |-)b<-|) 



(3) 



i.e. in the light of ^ tq = |+)(+| and ti = |-)(-|, where 
1+) = (|0) ± |1))/V2. Although manifestly correlated, state 
(O is fully classical since it is diagonal in the basis resulting 
from the tensor product between {\0}a, \1}a] and {\+}b, 
namely two local orthogonal bases lll3ll . Thus £)(po) - 0. 



Assume now that while A is well protected from the external 
environment so is not B, which is in contact with a dissipa- 
tive Markovian bath. In such a case, the system dynamics as a 
function of time t is fully described by an amplitude-damping 
channel This has an associated quantum map &p that 
transforms state po according to £p(po) - Eap(,E]^+ E\PqE^^, 

where = |0)b<0| + VT^|1)b<1| and Ei = V?I0>b<1| are 
the associated Kraus operators while p is a probability that 
grows with time t according to /:> = 1 - e^~^' {y is a relaxation 
rate). As anticipated, &p{po) evidently belongs to class (|2]i 
for any p (i.e. Vf) since £0 and £1 act on B only (hence en- 
tanglement never appears throughout). In Fig. 1(a), we plot 
D [cf. Eq. ([T]i] against p and the rescaled time jt as result- 
ing from numerical evaluation lll4ll . Dip), which is initially 
null as discussed above, at a first stage grows and eventually 
decays to zero as /:> — » 1, i.e. for f » y"'. This marks a pro- 
found difference between QD and entanglement in that a lo- 
cal non-unitary and memoiyless channel is able to create QCs 
previosly fully absent. Remarkably, the dissipative dynamics 
is merely detrimental to the classical correlations C - I -T> 
|3], where I -S{pa)+S{pb)-S{p) is the mutual information 
1 21]. C indeed exhibits a monotonic decay vanishing for p^l, 
i.e. yf^ 00, as shown in Fig. 1(b). 

In the following, we make this result rigorous by analyti- 
cally deriving Dip) and C{p) so as to reproduce Fig. 1. To 
this aim, we first explicitly derive in a compact form the QD 
of states (|2]i when tq and t\ have the same purity. Next, we 
present a picture in terms of trajectories in the Bloch sphere 
clearly highlighting the physical mechanism that causes QD 
to necessarily grow in the present process. 

Among states (O a prominent instance is the resource state 
for the B92 quantum cryptography protocol [2], which reads 
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FIG. 2. (Color online) (a) Function S{x,y) for |x| < I and |>i < 1. 
(b) and (c) Function 6{x, y) for (x, y) e D as given by Eq. l|9j when 
.So = = 1 and (p = 7t/3 (b), if = 2n/3 (c). In either case, the red region 
on the x-y plane represents the corresponding D. 



PB92 = 5 (|0)A<Oi® |0>B<0| + |1>A<1| ® |+>B<+|). Similarly to the 
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popular BB84 [2] the above is among those quantum protocols 
where the exploited quantum resource is not entanglement, 
which is is fully absent. Rather, it harnesses one-way QCs 
stemming from the non-distinguishability of states |0) and |+) 
I I3II . Later on, we will indeed show that pB92 possesses the 
maximum allowed QD within family ^ with tq and ti having 
equal purities. Also, states such as (|2]i can allow for quantum 
locking (see Iil5i1 and references therein). To our knowledge, 
the literature lacks in explicit formulas for the QD of (|2). We 
will thus carry out an ab initio calculation. 

To this aim, our first step is to express the single-qubit states 
To and Ti in Eq. ^ through the Bloch-sphere representation 
as Tj = (1 +S/ ■cr)/2, where 1 and cr = {0-1,0-2,0-3) are the 
usual identity and Pauli operators, respectively, while s, is the 
Bloch vector corresponding to r, (/ = 0, 1). Without loss of 
generality we can assume that io.v = ■sov - 0, ^O; = ■so and iiy = 0, 
i.e. in the Bloch sphere So and Si lie on the X-Z plane with 
So along the Z-axis (5, = |s,| < 1). Indeed, one can reduce 
the problem to such a case by applying a suitable single-qubit 
rotation, which cannot affect the QD like any local unitary 
operation Jst). Using this along with |0)(0| — (1 +0-3)72 and 
|1)(1| = (]1- 0-3)72, state (|2]) can be arranged as (from now on 
we drop subscripts A and B) 



p=i [1®1 + ]1( 



•(fllO-i +fl30-3) + 0-3®(/7iO-i -l-i'30-3)] (4) 



with 



ai- 



^isin^ 



(33 : 



.^Cff.S'iCOS ip 



So-S\ COS ifi 



(5) 



where we have carried out the replacements si^- - iisini^, 
= si cos (fi {(f is the angle between So and Si). Next, in the 
light of (HJ we need to calculate how p is transformed when 
a generic Von Neumann measurement is performed on sys- 
tem B, i.e. pk, as well as the associated probability pk- As 



3 



in III6II we use the property that the projector corresponding 
to any such measurement can be expressed as 1 (8> Bk with 
= VUuVK where = \k){k\ (k = 0, 1) and V is a generic 
one-qubit unitary. We thus expand p^.pj. -{i® Bt.)p{t® Bf^) as 

p^p,,^ {i®v){i®nk){t®v^)p{t®V){t®nk){i®v^). (6) 

The unitary V'' tranforms each PauH matrix cr, according to 
yV,V — Yj'p\ VijCj, where vij are real numbers satisfying the 
constraint L'j^.+i'2y+i'3^- = 1 , for any y = l,2,3. Also, for A: = 0, 1 
n^cririi =nj:cr2rij; =0 while 11^:0-311^: = fkUk, where fo-l and 
fi--l. This along with (01 yield ^ in the simpler form 



PkPk^(i®V) 



jU<:1+VaO-3 



where 



Vk^fkibiVu+biUii) . (7) 



As yni V' represents a pure state and each Pauli matrix is 
traceless it is immediately checked that pk -^kl'^- Therefore, 
we obtain that pk - (^+VklnkO'^)l'^y^ky) whose eigenvalues 
are {l + Vk/fik)/'2 (each two-fold degenerate). The quantity to 
minimize entering the last term of ([U thus reads 



(8) 



k=Q,l 



where h{x) - -xlog2 x - (1 - jc) log2(l - :c) is the binary 
Shannon entropy function and on the l.h.s. we have high- 
lighted the dependence on variables un and 1133. The iden- 
tity i;j3 -n/23 + 1^33 = 1 (see above) yields that {^13,1/33) must 



fulfill v?, + V. 



1- 



^23 



< 1, i.e. they belong to the unit cir- 



cle. To work out the infimum of 6 we accomplish the linear 
transformation x - aiVi-i + a-iV-iT, and y = biV\^+b-iV-},-i- This 
way, ([8]) now becomes a universal function of x and y, which 
we call S{x,y), i.e. it no longer depends on parameters {a,, bi] 
specifying states tq and ti [cf. Eq. Q]. 

In Fig. 2(a) we plot 5 against x and y. Note that its actual 
domain of definition is the square having side length 1/ V2 and 
vertices at points (±1,0) and (0,+l) [this is because of the log- 
arithms of \+fkyl{l+fkx), see Eqs. ©-(IH])]. Also, 5 is an even 
function of x {y) for any set value of y {x) as is also clear from 
its dependance on x and y through Hk and Vk [see Eqs. (|7|i-(l8]l]. 
It takes value 1 for x - y - Q and decreases as the distance 
from the origin grows. Its concavity is minimum along the 
X-axis, where the function is fully flat, and maximum along 
the y-axis [see the sail-like shape in Fig. 2(a)]. As for the do- 
main b within which ^ is to be minimized, i.e. the region 
in the new reference frame corresponding to the unit circle 
Uj3+L'33 < 1 (see above), it is immediately checked through the 
inverse transformation that this is an elliptic region given by 



b^[x,y: Ax-+2Bxy+Cy- < 1} , 



(9) 



A^(^^V , B^ , , C^(^^) . (10) 

As is evident from Fig. 2(a), owing to the concavity of 6 its 
infimum under constraint (|9) necessarily lies on the boundary 
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FIG. 3. (Color online) D (a) and C (b) vs. .v and ip. 

of D, i.e. the ellipse obtained from (|9|l by turning the inequal- 
ity into an identity. While such ellipse is centered at x-y-Q, 
its own axes in general do not coincide with those defining 
the reference frame. Yet, in the case that B - [cf. Eqs. ^ 
and ( fTOb l the ellipse is not rotated with respect to the x- and 
y-axis. This circumstance physically occurs when tq and ti 
have the same purity, i.e. so - s\ [see Eq. (fTOlil. Henceforth, 
we will focus on such a case, which is enough for the scopes 
of the present work and allows for a prompt analytical deriva- 
tion of the infimum of (O under constraint (|9). Hence, by 
setting So - si - s coefficients ( fTOl ) become A - Z? = 0, 
C - l/r^., where Tj. - s\ sin((^/2)| and r,, - s\ cos((/j/2)| are the el- 
lipse semi-axis lengths along the x- and y-axis, respectively. 
Now, because of the shape of S(x,y) [see Fig. 2(a)] when 
the main axes and the ellipse own ones are collinear S{x,y) 
takes its minimum at the points (0, ir,,). This is evident from 
Figs. 2(b) and (c), where we plot the restriction of 6 to region 
D in a paradigmatic case such that rv > r, (b) and one featuring 
< r,, (c). The minimum of (H)) within region (|9]l, i.e. the last 
term in Eq. thus reads 

min5(x,3;) = ^(0,r,) = S(0,s|cosf|) = /![i^^|^] , (11) 

where we used that /ik - 1 +fkX and Vk -fky [cf. Eqs. (111)-®]. 

As for S{pb), the trace over A of (|4| is obtained as pB — 
(l+flicri+fl3cr3)/2, whose eigenvalues with the help of Q are 
found as (l±s| cos (/j/2|)/2. Hence, 5'(pb) = /i[(1+s| cos (p/2\)/2]. 
Using again Eqs. (01) and the eigenvalues of p are calcu- 
lated as (l+s)/4, each being two-fold degenerate yielding that 
S{p) = 1-h/i[(1+s)/2]. Using these results along with (fTTT i in the 
light of ([T]), we find the QD of any state ^ such that tq and 
Ti have the same Bloch-vector length s in the compact form 

^^^^ p+..|cos.,/2| J^^^ p+.|.sinW2| J_^ (12) 



As for the classical correlations C-I- 
5(]1a/2) = 1 we find that 



C=l 



-„[ 



lj£|sin^/2| 1 



£), using that S (pa) - 



(13) 



In Fig. 3 we plot D (a) and C (b) vs. s and if as given by 
Eqs. ( fT2] i and ( fT3] l. For a given s, i.e. for fixed purity, D 
solely depends on \(p - n/2\ exhibiting its maximum value 
when if - 71 12. As \(f-nl2\ grows the QD progressively de- 
creases until it vanishes at 1^ = 0, n. On the other hand, for 
a given if, D increases with the purity at a rate that grows 
as If approaches nl2. Within the present class of states, the 
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FIG. 4. (Color online) Trajectories of £,j(|+>(+|) (blue solid line) and 
£;)(|-){-|) (red dashed) on the XX plane of the Bloch sphere. 



QD thus takes its maximum value ©max - 0.202 [about three 
times larger than the maximum value attained in Fig. 1(a)] for 
s- \ and ip-Tij2, which correspond to the B92 resource state 
Pb92 introduced previously. This feature is in accordance with 
Ref. ITtIi (where a different QCs' measure was used). As for 
C, Fig. 3(b) shows that this is maximum for s - \ and ip-n 
and decreases when either of such two parameters is reduced. 

We are now in a position to provide a comprehensive ex- 
planation for the effect in Fig. 1. Using the above discussed 
Kraus decomposition of £>p it is straightforwardly found that 
g p(|+) <+|) = (1 + Vl-po"! +po-i)l2. Hence, iox = -s\x = 
-y/1 -p, while - s\~ - p [note that this immediately shows 
that 50 = ^1=^ for any /?]. Using this, we end up with 

s{p)^ yjl+pip-l), ip{p)^n-2axcisn[pl ^Jl-p) . (14) 

By replacing flAi in Eqs. (fT2l i and ( fT3] l, the functions D{p) 
and Cip) so obtained reproduce the plots in Fig. 1. 

To shed light on the physical mechanism behind the effect, 
in Fig. 4 we plot on the XZ plane of the Bloch sphere the 
parabolic trajectories of the Bloch vectors corresponding to 
Sp(\+}{+\). Initially, the two vectors have unit length pointing 
towards opposite directions of the X-axis. Hence, i = 1 and 
ip^TT, which yields that £) = and C = 1 [cf. Eqs. ^ and (O 
and Fig. 3]. When p = 1, i.e. for yf » 1, either state is mapped 
onto |0)(0|, i.e. the North Pole, which gives that s = I, ^ = 
and thereby D = C - [Eq. ^ indeed shows that the state 
becomes evidently uncorrelated] . At the intermediate stage, 
map Sp acts in a way that the two Bloch vectors move along 
symmetrical parabolic trajectories (see Fig. 4) yielding that, 
in particular, <(p <n. On the other hand, the purity shows 
a dip with s taking its minimum s^m - V3/2 for /;> = 1/2 in 
accordance with (fT4l i. In the light of Eq. (fT2l i and Fig. 3(a), 
this necessarily brings about that D>0 since 0<ip<n and the 
Bloch vector length, although shrinking to some extent, keeps 
finite throughout. This clarifies why discord must necessarily 



be created in the course of the process. In essence, the local 
dissipative channel transforms and |-)(-| so as to re- 

duce their distinguishability, which unavoidably gives rise to 
QCs without entanglement Js, 13). Note that the process non- 
unitarity is crucial since distinguishability as measured by the 
scalar product is unaffected by any unitary Isl UHj. Also, note 
that based on the definition of classically correlated states lll3ll 
the above reasoning shows that QCs will be surely created re- 
gardless of the specific measure chosen to quantify them. 

In summary, in this Letter we asked whether QCs can de- 
velop as a result of local non-unitary dynamics, an unattain- 
able phenomenon with entanglement. After deriving the QD 
for the class of involved states, we analytically proved that this 
indeed can occur for two qubits, initially in a fully classical 
state, under a local memoryless amplitude-damping channel. 
Also, we showed that the mechanism behind the QCs' birth 
can be readily grasped in the Bloch-sphere picture. 

All these phenomena are arguably not restricted to qubits. 
The generalization to continuous-variable systems in a way 
that local bosonic Gaussian maps now play the role of the 
amplitude-damping channel is under ongoing investigations. 

We thank R. Fazio and M. Patemostro for comments, 
and acknowledge support from FIRB IDEAS through project 
RBID08B3FM. 

Note added. -Dming completion of this work we became 
aware of a related manuscript lll9ll . 
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